1
От начальных значений к граничным условиям
MATH007Lesson 11
00:00
Представьте разницу между выстрелом снаряда (где результат зависит от начального угла и скорости) и натяжением высокопрочной проволоки между двумя небоскрёбами. В первом случае вы задаёте начальные условия и смотрите, куда он попадёт; во втором случае проволока должна должна приземлиться в определённом окне на втором здании. Этот переход от «шагового» движения к «связанному» движению определяет смену задач с начальными значениями (IVP) на задачи с граничными условиями (BVP).

Определение задачи с граничными условиями

Стандартная задача с граничными условиями второго порядка включает дифференциальное уравнение, определённое на интервале $[a, b]$, где состояние системы фиксируется на обоих концах. Математически это выражается как:

$y^{\prime \prime}=f(x, y, y^{\prime}), \quad \text{ для } a \leq x \leq b$

с граничными условиями Дирихле:

$y(a)=\alpha \quad \text{ и } \quad y(b)=\beta$

Ключевое различие

В отличие от задач с начальными условиями, которые требуют $y(a)$ и $y'(a)$ в одной точке, задачи с граничными условиями указывают значение $y$ в точках $a$ и $b$. Мы больше не знаем «начальный наклон» $y'(a)$; вместо этого мы должны определить траекторию, которая «соединяет точки», при этом удовлетворяя управляющему уравнению во всей внутренней области.

Существование и единственность (Теорема 11.1)

Хотя теорема Пикара–Линдёфа гарантирует локальную единственность для задач с начальными условиями, задачи с граничными условиями определяются глобальным поведением. Даже простое линейное ОДУ может не иметь решений, иметь одно единственное решение или бесконечно много решений в зависимости от длины области $(b-a)$. Уникальное решение гарантировано, если:

  • $f$, $f_y$ и $f_{y'}$ непрерывны на области.
  • $f_y > 0$ (это действует как «восстанавливающая сила», гарантирующая, что решение не уходит в бесконечность).
  • $|f_{y'}|$ ограничено некоторой константой $M$.

Практическое применение: прогиб конструкции

Рассмотрим конструктивную балку длиной $l$, подверженную равномерной нагрузке $q$ и горизонтальному растягивающему усилию $S$. Прогиб $w(x)$ определяется уравнением:

$\frac{d^2 w}{d x^2}(x)=\frac{S}{E I} w(x)+\frac{q x}{2 E I}(x-l)$

С граничными условиями $w(0)=0$ и $w(l)=0$. Здесь концы балки закреплены, и нам нужно найти кривую $w(x)$, описывающую физическую форму балки под напряжением.

🎯 Основная численная философия
Переход к задачам с граничными условиями требует нового численного инструментария. Мы не можем просто интегрировать вперёд, потому что начальный наклон $y'(a)$ — это неизвестный «угол стрельбы», который необходимо настроить до тех пор, пока мы не достигнем цели $\beta$ в точке $x=b$.